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ABSTRACT 



This paper represents a first attempt of embedding the 
Barycentric Celestial Reference System, the fundamen- 
tal relativistic reference system to be used for the model- 
ing of Gaia observations, into some cosmological back- 
ground. The general Robertson-Walker metric is trans- 
formed into local coordinates where the cosmological ef- 
fects are represented as tidal potentials. A version of 
a cosmological BCRS is then suggested to lowest or- 
der. The effects of cosmological background on the mo- 
tion of the solar system are estimated and found to be 
completely negligible. The relation to the de Sitter and 
Schwarzschild-de Sitter solutions is discussed. 
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(1) tidal forces due to any particular external body (e.g., 
nearby stars or the Galaxy), 

(2) effects of cosmological background. 

The tidal forces can be easily taken into account using the 
same theoretical framework as that used to construct the 
BCRS (e.g. Soffel et al. (2003)). These effects were esti- 
mated to be negligible even for GAIA accuracy and, thus, 
the corresponding terms have been dropped in the BCRS 
metric tensor. To include effects from the cosmological 
background a new approach will be necessary. 



3. THE ROBERTSON-WALKER METRIC IN 
THE VICINITY OF AN OBSERVER 



1. INTRODUCTION 



The re lativistic model for GAIA published by iKlioneJ 
is formulated in the first post-Newtonian approxi- 
mation under the assumption that the solar system is iso- 
lated. The relativistic reference system where the posi- 
tions and other parameters of the sources are described 
is the Barycentric Celestial Reference System (BCRS) 
of th e International Astronomical Union ( So ffel et a/.L 
This model can be refined in two principal ways: 
(1) one can account for various sources of gravitational 
fields generated outside the Solar system (cosmology, mi- 
crolensing, and gravitational waves); (2) one can opti- 
mize the post-Newtonian formulas to allow a faster nu- 
merical implementation of the model. Several projects 
aimed at such refinements are underway. Here we dis- 
cuss the imbedding of the BCRS into some cosmological 
background. 



2. THE BCRS AND EXTERNAL GRAVITA- 
TIONAL FIELDS 



The BCRS is constructed under the assumption that the 
solar system is isolated. One neglects herewith two kinds 
of effects: 



Let us first consider a general Robertson-Walker met- 
ric in the neighborhood of an observer, the size of the 
neighborhood being about the size of the solar system. 
Thus, assuming the usual conditions for homogeneity 
and isotropy the cosmological metric can be shown to 
coincide with the well-known Robertson- Walker metric 
which in spherical coordinates [x'^ = ct, p, 9, tf) reads 



ds'^ = -c^dt'^ + a^{t) 
dp" 



+ p^de^ sm^Odip^ 1 , (1) 



1 - fcp2 



with k = —1, or +1 being the curvature parameter By 
a simple transformation of the radial coordinate 



, 1 . 

p = r IH — kr' 
I 4 



this metric can be transformed into isotropic form 

ds^ = gap dx" dx^ , 
500 = -1, 
9oi = 0, 



(2) 



(3) 



where {xf' = ct, x^), r — are Cartesian coordinates. 
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It is clear that the following coordinate transformations 

2s 



a{t)r 



a{t)r 



(4) 



(5) 



Note that As and Bs are functions of T in (|6|l. The above 
formulas ( I10> for these quantities should be understood 
in such a way that t should be formally replaced by T 
in the right-hand side of each equation. E.g. AiiT) = 

^© . A program for Mathematica written by the 

authors allows one to perform the matching procedure to 
arbitrary order of r. 



can be used to transform the metric Q into a form (this 
form will be called local metric below) where the non- 
Minkowskian terms appear only as powers of radial dis- 
tance from the origin 



Thus, we have constructed the general Robertson- Walker 
metric in local coordinates where the cosmological ef- 
fects manifest themselves only as tidal potentials. 



Goo = -l + E^^(^) (-) ' 



Goi — 0, 



Gab = 5ab{l + Yl^Bs{T) (6) 

where — cT,X^ are the local coordinates, R — \X^\, 
and As, Bs, Cs and Dg are unknown functions to be de- 
termined from the tensorial transformation law: 



5eA(i,x) 



dX"" 



G^,(T,X). (7) 



The matching procedure allows one to derive the un- 
known functions in a unique way. Let us introduce the 
following notations: 



1 , , 



(8) 
(9) 



Then one has 



^1 =P2, 
1 



A2 = - {pt ~ 2pIp2 -pI + HpI- P2)), 
Bi^-^ipl + f), 

B2^-^{pt~Aplp2-2fpi~3f), 

C2 = \pi (2P2-/), 

d^ = \pI 

D2^^pI (2P2-/), 



(10) 



It is clear already from the dimensionality consideration 
that the coefficients As, Bs, Cs, Ds are homogeneous 
polynomials of the differentiation operator d/dt acting on 
a{t) and / weighted as a term of degree two. The degrees 
of the polynomials Ag, Bs, Cs, Ds are 2s, 2s, 2s — 1 and 
2s. 



4. A PARTICULAR CASE: THE FLAT DE SIT- 
TER SOLUTION 



In the particular case of a flat de Sitter universe 
(flat empty universe with a A term) one has a{t) = 
exp(c y/A/St) and fc = in Q, so that the metric reads 

500 = ~1, 
got = 0, 



(11) 



= 6,j exp I c \/ t 



It is also well known that the flat de Sitter solution can be 
written in another form 



ds' 



-A dT^ + A-^ dp^ + {d0^ + sin^ 



1 



3 ' 

which under the transformation 

2p 



R 



i+yr~iA7 

R 



P = 



1+ j^Ai?2' 



(12) 

(13) 
(14) 



can be re-written in isotropic form 

ds^ = G„/5 dX" dX^, 



Goo — — 
Go. = 0, 

Gab = Sa 



1 



i-2^R' 

1 + A i?2 



(15) 



{l + i-^ARr 
It is easy to see that il5i coincides with ^ in the case 
a{t) = exp y/A/Stj and fc = 0. 



It is also known that the transformation (iLemaitrei 119251 
iRohertsonl (19281 iTolmanL 11931 



P 



cxp 



cT 



(16) 
(17) 
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can be used to transform il2\ into (^J and vice verse. 
Inverting ( I16t -(I17> one gets 



p — ar 



1--Aa^r^j, (18) 
(19) 



with a(t) = exp(^c^/X/3ty Combining ([T8t-(fT9t 

with il4\ one gets the transformation between the two 
isotropic versions of the flat de Sitter metric il 1\ and 



r = i- ^ a/^ log I 1 - ^ Aa 



R 



2c V A 
2ar 



1- iAa2 r2 



(20) 
(21) 



The coefficients of the series expansions of (I20> -( I21> co- 
incides with the corresponding coefficients ilQ\ in (0}- 
This gives an additional check of the general formu- 
las given in Section|3l 



5. THE COSMOLOGICAL BCRS METRIC TO 
LOWEST ORDER 



The cosmological BCRS metric can be considered as a 
perturbation of the corresponding cosmological metric. 
The gravitational field of the central body (or a system of 
bodies, e.g. of the solar system) should be embedded into 
the cosmological metric. Neglecting the interaction of 
the cosmological fluid (including the cosmological con- 
stant) with the solar system matter one gets the following 
simple version of the cosmological BCRS metric: 

500 ~ -1 + § ^(^' ""^ " ^ + ^ ^^^^^ I''''' 

4 

g,j « % (l + ^ ^(^' ^) + ^ ^i(^) I'^l') ' (22) 

where w and are the usual BCRS metric potentials. 
Here we neglected: (1) higher post-Newtonian terms 
0{c~^) in goQ and in goi, and 0{c~*) in gij due to post- 
post-Newtonian (and higher order) effects from the so- 
lar system matter; (2) higher-order cosmological terms 
C'(|x|'*) proportional to Ag and Bg with s > 2, (3) any 
terms induced by the interaction between the solar system 
matter and the cosmological fluid. For the case when the 
only source of energy-momentum tensor is the solar sys- 
tem matter and the field equations contain the cosmolog- 
ical constant A Soffel & Klioner (2003) have explicitly 
demonstrated how the metric i22\ can be derived from 
the field equations in the corresponding approximation. 



6. THE LOCAL EXPANSION HYPOTHESIS 



The apparently simple question whether the cosmological 
expansion happens also locally (say, whether a hydrogen 



atom or the Solar system also expand) is a very compli- 
cated one and presents still an unsolved problem. Starting 
from Einstein himself different auth ors got differen t an- 
swers using different arguments (see Bonnor (2000") for a 
review of recent progress). Certainly, the answer to this 
question crucially depends on our model of the matter in 
the universe. 

The cosmological BCRS metric suggested above implies 
that the cosmological expansion has certain influence on 
the properties of space-time within the solar system. This 
is certainly true for the terms coming from the cosmolog- 
ical constant A or vacuum energy since this energy source 
is present everywhere. Note that recent cosmological ob- 
servations suggest that about 73% of the energy in the 
universe comes from that source. The applicability of 
the suggested BCRS metric to the 4% coming from the 
luminous matter and the 23% of the dark matter has to 
be investigated further. Ideas of cosmological averaging 
might be employed here (.Futamase. . 1 996i) . 



7. THE SCHWARZSCHILD-DE SITTER SOLU- 
TION 



The well-known Schwarzschild-de Sitter solution pro- 
vides additional arguments in justifying the cosmological 
BCRS metric. The standard form of the Schwarzschild- 
de Sitter metric resembles il2\ : 



ds^ = -A dT^ + A-^ dp^ + ( 

A = l -Kp\ 

P 3 



sm 



(23) 



This solution corresponds to the Schwarzschild solu- 
tion with mass m for the field equations with cosmo- 
logical constant. For A = this solution coincides 
with the Schwarzschild solution in standard coordinates, 
and for m = thi s agr ees with the de Sitter solu- 
tion M2\ . Robertson ( 1928) has shown that the metric 
( I23> under the transformations (in this section a{t) = 
exp(cA/A73i)) 



Sm^ /A / 1 ^ 2 am r 

T = t+ W-F' ^' 



F{x) 



c V 3 \ 14 

dx 



2 a{t) r 



(l-2;2)(a;2(l-a;2)2-|m2 A) 



p = a{t)r 1 + 



2a{t)r 
can be put into form 

ds^ — gap dx" dx^ , 



(24) 



,(25) 



(26) 



1 - 



500 



5oi = 0, 



2 a(t) r 



1 



2 a(t) r ^ 



<5y a{tf ( 1 + 



2a{t)r 



(27) 
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For A = this metric coincides with the Schwarzschild 
solution in isotropic coordinates, and for m = this met- 
ric gives the de Sitter solution in the form (II U . It is im- 
portant to note that Eqs. ( I24t -(I26> coincide with (II 8> - 
([19} for m 0. 

Metric (123 > written in isotropic or harmonic gauge rep- 
resents an exact solution of the cosmological BCRS in 
local coordinates (cosmological terms as tidal potentials 
only) for the particular case of the de Sitter metric with a 
spherically symmetric perturbation. This solution can be 
used to verify the cosmological BCRS metric in higher 
approximations. The details will be published elsewhere. 



8. DYNAMICAL COSMOLOGICAL EFFECTS IN 
THE SOLAR SYSTEM 



The cosmological terms in the suggested BCRS metric 
lead to an additional central force described by the dis- 
turbing function 

R^^Air\ (28) 
The parameter Ai can be calculated as 

Ai^P2 = - = -qH^, (29) 
a 

where H and q are the Hubble constant and the deceler- 
ation parameter, respectively. Using the current best esti- 
mates of the parameters q « —0.6 and H ^ 71 km/s/Mpc 
one gets 

« 3.2 X 10"^^ s"^ (30) 

This is the value of Ai at the present moment. The time 
dependence of this quantity is different for different cos- 
mological models. In the flat de Sitter universe Ai is time 
independent: Ai — ^ A — const. In general Ai, how- 
ever, is time dependent. At the distance of Pluto (r = 40 
AU) the maximal value of the disturbing cosmological 
acceleration amounts to Ai x 40 AU w 2 x 10^^"^ m/s^. 
Let us check how the cosmological perturbations under 
study change the orbit of planets. Since the disturbing 
function (I28> gives a central force, the osculating inclina- 
tion i, osculating longitude of the node fl, and osculating 
semi-latus rectum p — a{l — e^), a and e being the semi- 
major axis and the eccentricity, remain constant: 

i = const, (31) 
ft = const, (32) 
p = a (1 — e^) = const. (33) 

Since the time scale over which Ai is changing is much 
larger than typical orbital periods in the solar system, 
we might use an adiabatic approximation and put Ai = 
const for the calculations of the other osculating ele- 
ments. In this case one gets that the secular changes of 
averaged argument of perihelion oJ and averaged mean 
anomaly Mq read 



AU = Ai rad/revolution, (34) 

7 + Sp^ 

AMo = ^— Ai p2 rad/revolution. (35) 
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Here and below P is the period of motion of the body 
under consideration. There are no secular effects in the 
semi-major axis and eccentricity and the amplitudes of 
periodic effects in these elements may amount to 

Ae = i— (36) 

4 TT^ 

Aa/a=-^A, = ^A,P\ (37) 

This shows that the secular perturbations in the argument 
of perihelion uj and in the mean anomaly M for Pluto are 
about 10^^ microarcsecond per century and three orders 
of magnitude smaller for Mercury. The amplitudes of the 
periodic effects in a and e are about 10^^^ for Pluto and 
2Q-23 Jqj. Mercury. Therefore, these cosmological per- 
turbations are completely negligible for the solar system 
dynamics over times much less than the age of the Uni- 
verse. 



9. FURTHER WORK 



Let us mention two almost independent problems which 
should still be solved: 

(1) A better version of the cosmological BCRS metric 
should be constructed where the non-linearity of the field 
equations is taken into account. 

(2) Parameters (e.g., parallaxes of distant sources) ob- 
tained while processing observations with the standard 
(non-cosmological) BCRS should be related to the pa- 
rameters meaningful in cosmological context. Let us note 
here, that the parallax distance defined as dp ~ a{t) r k. 
R agrees with the coordinate distance in our local coor- 
dinates. 
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